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Abstract 

We propose to undertake a model-independent analysis of the inclusive decay rates and 
distributions in the processes B — ► X s ~f and B — > X s £ + £~ (B = i? ± or We show how 

measurements of the decay rates and distributions in these processes would allow us to extract 
the magnitude and sign of the dominant Wilson coefficients of the magnetic moment operator 
TObSiCT MI/ 6_RF A " / and the four-fermion operators (si^&l)^^) and (sz,7 jU &l)(^7 m 7 5 ^)- 



1. The Decay B 
Experiment 



X s j in SM and 



The measurements of the decay mode B — > K*j, 
reported last year by the CLEO collaboration |Q, 
having a branching ratio B{B — > K*-f) = (4.5 ± 
1.0 ± 0.9) x 10~ 5 , and the inclusive decay B — » X s j, 
reported at this conference [^| with a branching ratio 
B{B -> X sl ) = (2.32 ± 0.51 ± 0.32 ± 0.20) x 10~ 4 , 
have put the physics of the electromagnetic penguins 
on an experimental footing. In the standard model 
(SM), these transitions are dominated by the short- 
distance contributions and provide valuable information 
about the top quark mass and the Cabibbo-Kobayashi- 
Maskawa (CKM) weak mixing matrix elements Vt s Vtb 
||. The rapport between the SM and experiment may 
be quantified in terms of the CKM matrix element ratio 



0.62 < 



Vts_ 

V cb 



< 1.1 



(1) 



which is consistent with unity, resulting from the 
unitarity constraints. Alternatively, one can set 
Vts/V c b — 1 to obtain from the CLEO measurement 
bounds on the Wilson coefficient G\(mi,) of the effective 
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magnetic moment operator. Using B(B — > A s 7) from 
one obtains 

0.22 < |C 7 (to 6 )| < 0.30. (2) 

Using, however, the 90%-confidence-level range from the 
CLEO measurement B(B -> X s j) = (2.31 ± 1.1) x 10~ 4 
and the theoretical calculation for B(B — > X s j) from [Q 
we obtain 

0.19 < \C 7 (m b )\ < 0.32. (3) 

We also remark that the photon energy and hadron mass 
spectra measured by CLEO are in good agreement with 
the SM-based calculations in ||. The bound (§) can be 
used to constrain the non-SM contribution to the decay 
rate B(B — > X s j) as discussed in these proceedings 



2. Motivation for a Model Independent 
Analysis of Rare B Decays 

The determination of |CV(TOfc)| from the inclusive 
branching ratio B(B — » X s ~/) is a prototype of the kind 
of analysis that we would like to propose here to be 
carried out for the rare B decays in general and for 
the semileptonic decays B — > X s £ + i~, in particular. 
First steps towards a model-independent analysis of the 
FCNC electroweak rare B decays involving these decay 
modes have recently been proposed in to which we 
refer for details and references to other related work. 
Here, we summarize the main assumptions and results. 
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The main interest in rare B decays is to measure 
the effective FCNC vertices in order to test the 
SM and search for new physics. We have argued 
that with some plausible assumptions these vertices 
can be parametrized through a limited number of 
effective parameters, which govern the rates and shapes 
(differential distributions) in rare B decays B — > X s -f, 
B -> X s ^~and B s -> £+£-. The search for physics 
beyond the SM in these decays can be carried out 
in terms of three effective parameters, CV(/z), Cg(fi) 
and Cio(/i), characterizing the strength of the magnetic 
moment and two four-fermion operators (sl7^&l)(^7' 1 ^) 
and (sl7 a ,6l)(^7 m 7 5 ^). This can then be interpreted 
in a large class of models. The presence of non-SM 
physics may manifest itself by distorting the differential 
distributions in B — > X s £ + l" . Some possible examples 
of such distortions have been worked out in § . Here we 
present profiles of the Wilson coefficients in the best- 
motivated extensions of the SM, namely the Minimal 
Supersymmetric Standard Model (MSSM). 

Our analysis is based on an effective Hamiltonian of 
the form 



n eff (b 



sX) = " vf At 



10 

E 

i=l 



Ci(n)Oi(n) . (4) 



where X stands for qq, 7, gluon and £ + £~ and At = 
V t *V t b- The operator basis O1...10 is given in || and is 
the same as in the SM, thereby restricting our analysis to 
cases, in which the effective Hamiltonian may be written 
as (|). 

3. Analysis of the Decays B — ► X s j and 
B -> X s £+£~ 

The experimental quantities we consider in this paper 
are the following: (i) Inclusive radiative rare decay 
branching ratio B(B — ► X s ~f); (ii) Invariant 
dilepton mass distributions in B — ► X s £+£~; 
(iii) Forward- backward (FB) charge asymmetry 
A(s) in B -> X s £+£-. 

The FB asymmetry A{s) is defined with respect 
to the angular variable z = cos 9, where 8 is the 
angle of the £ + with respect to the 6-quark direction 
in the centre-of-mass system of the dilepton pair. 
It is obtained by integrating the doubly differential 
distribution d 2 B/{dzds)(B -> X s £+£-) g: 



A{s) 



dz 



d 2 B 

dz ds 



dz 



d 2 B 

dz ds ' 



(5) 



where s — (pi + P2) 2 /m 2 and p\ and pi denote, 
respectively, the momenta of the £ + and £~ . 

We remark that the decay rate B(B —* X s j) puts 
a bound on the absolute value of the coefficient CV(/x). 



However, the radiative B decay rate by itself is not able 
to distinguish between the solutions CV(/i) > (holding 
in the SM) and the solutions CV(/u) < 0, which, for 
example, are also allowed in the MSSM as one scans 
over the allowed parameter space. We recall that the 
invariant dilepton mass distribution and the forward- 
backward asymmetry in B — > X s £ + £~ are sensitive 
to the sign and magnitude of Cj(p) ||. Using H e ff 
given in (^) , one obtains for the dilepton invariant mass 
distribution 



Ay 



w(s) 



— = B — 

ds s Un 2 \V cb \ f(m c /m b ) 



(G) 



(\C 9 + Y(s)\ 2 + C 2 10 ) ai(J> s ) 
4 

--C 2 a 2 (s,m s ) + 12a 3 {s,m s )C 7 {C 9 + Re Y(s)) 
s 



where the auxiliary functions a, depend only on the 
kinematic variables and Y(s) depends on the coefficients 
G\ , • • • , Cg of the four quark operators (see (8) ) . 

The corresponding differential asymmetry as defined 
in (H) is 



A(s) = -B. 



3a 2 



1 



:W 2 (s)C w 



8tt 2 f(m c /mb) 
[s(C 9 + ReY(s))+4C 7 (l 



(7) 



rh 2 



)] 



We first present the partial branching ratio B(As) 
and partial FB asymmetry A(As), where As defines 
an interval in the dilepton invariant mass. In order 
to minimize long-distance effects we shall consider 
the kinematic regime for s below the J/tp mass (low 
invariant mass) and for s above the mass of the tp' (high 
invariant mass). Integrating ^ over these regions for 
the invariant mass one finds 

B(As) = A(As) (Cl + C 2 ) + B{As)C 9 + C(As), (8) 

where A, B and C are fixed in terms of the Wilson 
coefficients C\ ■ ■ ■ Cq and C7. For the numerical analysis 
we use nib = 4.7 GeV, m c = 1.5 GeV, m s = 0.5 GeV. 
The resulting coefficients A, B, and C are listed in Table 
[l] for the decays B — > X s e + e~ and B — > 

For a measured branching fraction B(As), one 
can solve the above equation for B(As), obtaining a 
circle in the C9-C10 plane, with centre lying at C 9 = 
B{As)/{2A(As)) and C 10 = 0. The radi us of this 
circle is proportional to y / B(As) — B m i n (As), where the 
minimum branching fraction 



B min {As) = C{As) - 



B 2 {As) 
AA(As) 



(9) 



is determined mainly by the present data on B — > X s "/, 
i.e. by |C 7 |. 
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As 


cv 


A(As) 


B(As) 


C(As) 
£ = e 


C(As) 
£ = n 


4mj < s < 

Am 2 < s < m 2 J/qp 
m 2 <s< (l-m?) 
mr^i < s < (1 — m 2 ) 


+0.3 
-0.3 
+0.3 


2.86 
2.86 
0.224 


-5.76 
-20.8 
-0.715 


84.1 
124 
0.654 


76.6 
116 
0.654 


-0.3 


0.224 


-1.34 


2.32 


2.32 



Table 1. Values for the coefficients A(As), B(As) and C(As) 
(in units of 10~ 8 ) for the decay B -+ X s £+£-. 



As 


C 7 


a(As) 


/9(A«) 


4m ?< s < m j/* 
4mf < s < m 2 /4i 

m 2 < s < (1 - mf) 
"V < s < C 1 - m ?) 


+0.3 
-0.3 
+0.3 
-0.3 


-6.08 
-6.08 
-0.391 
-0.391 


-24.0 
55.4 

0.276 
1.37 



Table 2. Values for the coefficients a(As)and /3(As) ( in units 
oflO- 9 ). 



To further pin down the Wilson coefficients, one 
could perform a measurement of the forward-backward 
asymmetry A. Integrating (||) over a range (As) yields 



A(As) = do (a(As)C 9 + /3(As)) 



(10) 



For a hxed value of A(As), one obtains hyperbolic 
curves in the C9-C10 plane; like the coefficients A, B and 
C, the parameters a and (3 are given in terms of the 
Wilson coefficients C\---C§, Cj and As; their values 
are presented in Table 

Given the two experimental inputs, the branching 
fraction B(As) and the corresponding asymmetry 
„4(As), one obtains a fourth-order equation for the 
Wilson coefficients Cg and Cio, which admits in general 
four solutions, which can be plotted as contours for 
a fixed value for the branching fraction B(As) and 
the FB asymmetry _4(As). The possible solutions 
for Cg and Cio are given by the intersections of the 
circle corresponding to the measured branching fraction 
and the hyperbola, corresponding to the measured 
asymmetry. Details can be seen in || . 

We stress that the spectrum itself is very sensitive 
to the values of the Wilson coefficients and to the sign 
of Cj. In fig. [l]we plot the various contributions to the 
spectrum, for positive and for negative C7. 

In a similar way, it may become possible to 
measure also the differential asymmetry (^|) ; the various 
contributions to A(s) are shown in fig. |2| 

4. Model Predictions for the Wilson 
Coefficients 

As an illustrative example we shall consider here the 
MSSM; we shall show, how the predictions for the 
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Figure 1. The dependence of the invariant- mass spectrum on 
the Wilson coefficients. Solid line: SM. Long-dashed line: 
CV — > — CV, with other coefficients retaining their SM values. 
Short-dashed line: The contribution of Cio onrv - Dotted line: 
Cio = 0, with other coefficients retaining their SM values. 
Dash-dotted line: same as for the dotted one, but with 
C7 = —0.3. The vertical lines indicate the location of the J/\P 
and W' resonances. 
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Figure 2. The dependence of the differential FB asymmetry on 

the Wilson coefficients. Solid line: SM. Long-dashed line: 

Cr — > — CV, with the other parameters retaining their SM values. 

The vertical lines indicate the location of the J/ty and \t' 

resonances. 



Wilson coefficients C7, Cg, and Cio are altered in this 
model. 

Once the gluino contributions (as well as the 
analogous ones from neutralino exchange) are neglected, 
the flavour violation in the supersymmetric models 
is completely specified by the familiar CKM matrix. 
The one-loop supersymmetric corrections to the Wilson 
coefficients C7, Cg, and Cio are given by two classes 
of diagrams: charged-Higgs exchange and chargino 
exchange |To|| . 

The charged-Higgs contribution is specified by two 
input parameters: the charged-Higgs mass (tojj+) and 
the ratio of Higgs vacuum expectation values (i»2 / v i = 
tan (3) . This contribution alone corresponds to the two- 
Higgs doublet model which has also been considered in 
§■ 

In addition to the diagrams with charged-Higgs 
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exchange, the MSSM leads also to chargino-mediated 
diagrams. The chargino contribution is specified by six 
parameters. Three of them enter the 2x2 chargino mass 
matrix: 

( M m w V2sinp\ 

Following standard notations, we call tan/3 the ratio 
of vacuum expectation values, the same that appears 
also in the charged-Higgs sector, and M, \i the gaugino 
and higgsino mass parameters, subject to the constraint 
that the lightest chargino mass satisfies the LEP bound, 
m+ > 45 GeV. The squark masses 

m? 2 = fh 2 + m 2 ± Amm (12) 

contain two additional free parameters besides the 
known mass of the corresponding quark m q : a common 
supersymmetry-breaking mass fh and the coefficient A. 
The last parameter included in our analysis is a common 
mass mr for sleptons, all taken to be degenerate in mass, 
with the constraint mj > 45 GeV. Therefore the version 
of the MSSM we are considering is defined in terms of 
seven free parameters. 

We have computed the Wilson coefficients in 
the MSSM and then varied the seven above-defined 
parameters in the experimentally allowed region. The 
results of our analysis are presented in fig. || which 
shows the regions of the C9-C10 plane allowed by 
possible choices of the MSSM parameters. The upper 
plot of fig. |3] corresponds to parameters which give rise to 
positive (same sign as in the SM) values of CV, consistent 
with experimental results on b — ► 57 (0.19 < CV < 0.32), 
while the lower plot corresponds to values of C7 with 
opposite sign (—0.32 < C7 < —0.19). We also show 
how our results are affected by an improvement in the 
experimental limits on supersymmetric particle masses, 
as can be expected from the Tevatron and LEP 200. 
Fig. H also shows the Cg — C10 regions allowed by the 
MSSM if the further constraints m H + > 150 GeV, 
mj , m x + , mj > 100 GeV are imposed. 

The regions shown in fig. || illustrate the typical 
trend of the supersymmetric corrections. If supersym- 
metric particles exist at low energies, we can expect 
larger values of C10 and smaller (negative) values of 
Cg than those predicted by the SM. This is the general 
feature, although the exact boundaries of the allowed 
regions depend on the particular model-dependent as- 
sumptions one prefers to use. However, the most inter- 
esting feature of supersymmetry is that solutions with 
negative values of C7 are possible and are still consis- 
tent with present data. Moreover, values of the other 
two coefficients Cg and C10 sufficiently different from 
the SM are allowed, leading to measurable differences 
in the decay rates and distributions of B — > X s £ + £~ 
and B s -► £+£-. 
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Figure 3. The region in the C9-C10 plane obtained by varying 
the MSSM parameters. The upper (lower) plot corresponds to 
solutions that satisfy the b — > sj experimental constraint with 
positive (negative) C7 given in eq. (El) and the present bounds 
{m H + > 80 GeV, rht , m x + , > 45 GeV). The smaller areas 
limited by the short-dashed line correspond to the region of the 
MSSM parameter space that will survive an unsuccessful search 
for supersymmetry at the Tevatron and LEP 200 {m H + > 150 
GeV, m h m x+ , m; > 100 GeV). 
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